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ABSTRACT 



We explore the initial conditions for fragments in the extended regions (r > 50 
AU) of gravitationally unstable disks. We combine analytic estimates for the 
fragmentation of spiral arms with 3D SPH simulations to show that initial 
fragment masses are in the gas giant regime. These initial fragments will have 
substantial angular momentum, and should form disks with radii of a few AU. 
We show that clumps will survive for multiple orbits before they undergo a 
second, rapid collapse due to H2 dissociation and that it is possible to destroy 
bound clumps by transporting them into the inner disk. The consequences 
of disrupted clumps for planet formation, dust processing, and disk evolution 
are discussed. We argue that it is possible to produce Earth-mass cores in the 
outer disk during the earliest phases of disk evolution. 



Keywords: DISKS; PLANETARY FORMATION; PLANETS, MIGRATION; 
PLANETESIMALS 
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1 Introduction 



Protoplanetary disks are likely to be massive during their initial phases of evo- 
lution. Collapse calculations demonstrate that gravitationally unstable disks 



do form (e.g., Vorobyov and Basu, 2009) and that their evolution includes 



phases of strong gravitational instability (e.g., Vorobyov and Basu, 2006). In 
the inner disk (r < 50 AU), cooling times relative to local dynamical times are 
too long for the instability to result in fragmentation, and the disk reaches a 



self-regulating state (e.g., Boley et al. 2006 Durisen et al. , 2007). There may 



be some exceptions due to changes in disk chemistry (Mayer et al. , 2007), but 



these situations require further study. In contrast, fragmentation in the outer. 



extended disk (r > 50 AU) becomes quite possible if the Toomre (1964) Q can 



be driven toward unity by, e.g., mass loading (Boley, 2009). Whether these 



fragments typically produce brown dwarfs (Stamatellos et al. , 2007 Stamatel 



los and Whitworth, 2009) or gas giant planets is a topic of debate. However, 



as we argue here, clumps that become destroyed can be just as important to 
planet formation and disk evolution as clumps that remain bound. 



In this paper, we study disk fragmentation conditions and make estimates for 
initial fragment masses; we describe the initial angular momenta of fragments; 
and we discuss the consequences of clump disruption for planet formation. We 
present simulation data and describe toy models in section 2. In section 3, 
we derive proper estimates for fragment masses, and find consistency between 
our estimates and simulation data. We estimate the expected initial angular 
momenta of clumps in section 4, which are also in rough agreement with 
simulation data. We use a polytropic model in section 5 to estimate the time 
during which a clump could be tidally disrupted, including accretion effects. 
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These sections culminate to show that clump disruption is a real possibility 
during disk evolution. For section 6, we speculate on the consequences of 
disrupted clumps for dust processing, core formation, and the FU Orionis 
phenomenon. Our conclusions are given in section 7. 



2 Models 



We present two parameterized models. Model A is used to show what we 
expect for fragmentation around A stars. We assume a temperature profile 

a mean molecular weight fi = 2.3, and the epicyclic 
frequency k ^ Q ^ f^Kepierian for central star M = 1.5Mq. The surface density 



S is given by the Toomre (1964) parameter Q = Csk/{'kGT,), where Cg is 
the local sound speed. The entire disk does not need to be gravitationally 



unstable {Q < 1.7, Durisen et al. , 2007), and the total disk mass for the 
model could vary considerably depending on the size of the low-Q region. If 
Model A's disk has a Q = 1.5 between r ~ 100 and 150 AU, the mass in 
this region would be about 0.15 Mq. Our second model. Model M, represents 
a disk orbiting an M star, embedded in an envelope with a 30 K irradiation 
temperature. For Q = 1.5, the mass contained between r ~ 100 and 150 
AU is 0.04 M0. In addition to these parameterized models, we present a 
simulation with initial conditions (ICs) based on the simualtion "SIMA" from 



Boley (2009) just before the disk fragments. We refer to this model for the rest 



of the paper as SPHSIM to avoid confusion with SIMA and the analytic models 
described above. In SPHSIM, a O.3M0 star is surrounded by an r ~ 400 AU 
disk that is accreting mass from an envisaged envelope at ~ IO^^Mq yr^^. This 
accretion rate is consistent with what one expects in a protostellar collapse 
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with a background temperature of 30 K. At the time our simulation begins, the 
disk has a mass of ~ O.IOM©. The inner radius is at r ~ 18 AU and an outer 
radius of r ~ 510 AU. We do not include additional mass loading, so SPHSIM 
represents a last-burst scenario, i.e., the final phase of fragmentation that this 
disk is likely to experience. The data were interpolated from the CHYMERA 



(Boley, 2007) cylindrical grid to an SPH realization via a density- weighted 



Monte Carlo sampling. 



SPHSIM is run using GASOLINE (Wadsley et al. 2004), which is a multipur- 



pose code designed to model structure at various scales, e.g., planetary disks 
as well as cosmological structure formation. One million particles are followed, 
with a particle mass ~ 2 x 10^'^ Mj. The spline force softening is set to 2 AU 
for the star and 0.5 AU for all other particles. The radiative cooling algorithm 
is the same as that described in [Boley (2009), but adapted to SPH. The cool- 



ing is calculated from V ■ F = — (367r)-'^/^s^^cr(T'^ — T-^J(At + 1/Ar)^^, where 
s = (m/pY^^ and At = sup for the local opacity k, particle mass m, and den- 
sity p. The factor SGtt comes from defining the radius of a resolution element 
as r = (3m/ {ATxp)Y^^ and setting the ratio between the radiative flux and the 
divergence of the fux to be the area over the volume of the resolution element. 



D'Alessio et al. (2001) opacities are used, with a 1 /xm maximum grain size. 
The irradiation temperature Ti„ = 30 K everywhere. This cooling approxi- 
mation is good for the outer disk regions, where midplane optical depths are 
< 1. The disadvantage to this approximation is that it neglects the effects of 
radiation from a collapsing clump on its surrounding medium. 

In addition to the radiative cooling algorithm described above, GASOLINE 



has been augmented with the Read et al. (2009) OSPH modifications in or- 



der to address the SPH limitations outhned, most recently, by [Agertz et al. 
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(2007). We find that OSPH is superior to other suggested modifications (e.g., 



Price, 2008) because (1) OSPH works for a self-gravitating fluid, (2) recovers 
the correct timescale for the growth of the Kelvin-Helmholtz instability with- 
out introducing new free parameters, and (3) allows for convergence testing 
of the hydrodynamics by increasing the number of neighbors without being 
compromised by the tensile instability. We also note that GASOLINE uses a 
fixed gravitational softening, while the SPH smoothing length h is variable, 
defined by 32 neighbors. Using a variable softening length can induce fluctu- 
ations in the potential energy of particles, which inevitably leads to errors in 
energy conservation. On the other hand. Nelson (2006) showed that fixed soft- 
ening could lead either to either enhanced or suppressed clumping. Although 
there are ways of improving energy conservation that would permit the use of 



variable softening lengths ( Price and Monaghan , 2007 ) , this has yet to be im 



plemented in GASOLINE and represents a future code development project. 
Instead, we have chosen to ensure that energy is conserved and set the mass 
resolution such that the gravitational softening length is larger than the typi- 
cal SPH smoothing length in dense structure, e.g., spiral arms (see discussion 



m 



Mayer et al. , 2004). The fragments that form in SPHSIM (below) have a 
median softening to smoothing length ratio of about 4, so SPHSIM should 
err on the side of slowing the collapse of fragments. Moreover, GASOLINE 
showed highly satisfactory results in the Wengen Test 4 comparison project, 
where a highly unstable disk was followed in detail by several codes. In particu- 
lar, Gawryszczak & Mayer (2008) reported strikingly good agreement between 
GASOLINE and the FLASH AMR code when following fragmentation in a 
self-gravitating disk. 

Three clumps form in SPHSIM (CI, C2, and C3), and their initial properties 
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are listed in Table 1. We computed the number of particles per Jeans mass 
according to Mjeans/Mparticie = 2.92c^/(MparticieG^/Vj,{L) ) souud speed Cs 
and peak density Ppeak in the fragmenting spiral arm. Fragmentation is nu- 



merically well- resolved, with Mjeans/A^particic ~ 2000, satisfying the Bate and 



Burkert (1997) criterion (see also |Nelson[ |2006D . Particles are flagged as frag- 



ment members using SKID (Stadel, 2001), which groups particles according 
to density gradients and then completes an iterative unbind for each particle 
in the group. When comparing the results to SIMA, it should be noted that 
Boley (2009) reported a clump mass of 20 Mj for the end of the simulation. 
This estimate included mass growth for about 1.75 orbits, so it does not rep- 
resent the clump's initial mass, which is between 4 and 5 Mj. Boley expressed 
reservation in the Letter about accepting the final masses from his simulations 
because, e.g., the radiative effects of the clump on its surroundings were not 
modeled, and the resolution was too low to follow the evolution of the clump 
itself. In addition, SIMA only forms one clump, while SPHSIM forms three. 
Although integrated or time-averaged quantities between two realizations of 
a simulation should give comparable answers, detailed structure, especially 
clumps, are extremely sensitive to initial conditions (see Wengen 4 Compar- 
ison ProjectQ. The difference in the number of fragments between SPHSIM 
and SIMA are not considered by us to be failures of either model. 



3 Clump Mass 

In this section, we calculate initial clump masses for unstable disks. Why do 
we care about initial clump masses? Although clumps are likely to accrete, the 

1 www.astrosim.net 
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evolution of the system will depend on its initial state. For example, a non- 
rotating, one Jupiter-mass clump will contract for about a fewxlO^ yr, until 



molecular hydrogen dissociates, causing rapid collapse (e.g., Tohline, 2002). 
We refer to this contraction timescale as the clump's primary contraction 
time (ti). For five Jupiter masses, Ti is reduced to a few xlO'' yr, and for 
an initial mass of ten Jupiter masses, the clump will collapse in < 10^ yr 



(Helled and Bodenheimer, 2009). A factor of ten in the initial mass can change 
the clump's initial evolution timescale by a factor of about 100. The spatial 
scale of fragmentation, which is related to the initial mass, also determines 
the rotational angular momentum (next section). Once the core of the clump 
collapses, this angular momentum barrier should lead to the formation of a 
circumplanetary/brown dwarf disk. This disk is expected to control the long- 
term accretion history of the clump, as material entering the Hill sphere of 
the collapsed planet will have angular momentum from disk shear. In order 
to avoid confusion between fragmentation, i.e., the formation of the clump, 
and the rapid collapse that follows dissociation of H2, we refer to the latter as 
dissociative collapse. 



From a practical standpoint, initial masses can be well-constrained using 
global simulations, while the subsequent evolution can only be modeled poorly 
at this time. Simulating clump growth requires resolving convection, the pho- 
tosphere, chemical changes, the clump's effect on the surrounding medium, 
the core/disk transition region, and the gas flow into the Hill sphere. Even 
at our resolution of about 5000 particles per Jupiter mass, this is a daunting 
task, and best addressed by high-resolution simulations of individual clumps. 
For all of these reasons, we argue that constraining initial clump masses is fun- 
damental to understanding the fragmentation process in disks and gas giant 
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planet formation. 

What do we expect for initial clump masses? In order to calculate the fragment 
mass, we need to know the local surface density and size scale for fragmenta- 
tion. A back-of-the-envelope estimate is the Toomre mass 

Mr = 7r(AT/2)^S (1) 



e.g., Nelson, 2006), where the Toomre wavelength. 



Xt = 2c^/(GE), 



(2) 



is roughly the most unstable radial wavelength for local sound speed Cg and 
smooth surface density S. The surface density can be calculated from the 
Toomre parameter using S = Cs^/{ttG (Q)), and this value can be used to 
find the Toomre mass for a given Cg and (Q), where we have taken k, ^ Q. 
The brackets are used to denote that {Q) is a smooth, axisymmetric quantity. 
We can rewrite the Toomre wavelength 

XT = 2n{Q)fJ,H, (3) 
where the local scale height 

H = cJUefg^). (4) 



and /e and fg are shape factors of order unity that depend on the equation 
of state and on self-gravity effects, respectively. Setting {Q) = 1 and calcu- 
lating a Toomre mass based on the most unstable radial wavelength for the 
unperturbed axisymmetric disk (Ay) includes material over a radial extent 
that is 2ttH. Gravitational instabilities can produce strong spiral waves and 
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local density perturbations, making such use of the axisymmetric measure At 
over large radial scales in the non-linear regime, where fragmentation occurs, 
extremely dubious. The Toomre mass is strictly an estimate for the disk mass 
that becomes incorporated into one wavelength of the resulting spiral waves. 
Fragmentation is best described in the context of spiral arms. 



Instead of using measurements that correspond to an axisymmetric disk, we 
use length scales and surface density perturbations that are appropriate for 
spiral arms. The radial extent of fragmentation can be estimated using the 



results of Durisen et al. (2008, hereafter DHP2008), who used the virial the- 
orem to show that a disk, under isothermal conditions, is most susceptible 
to fragmentation within a region 5r from the corotation of a spiral wave. 
They demonstrated this behavior using isothermal hydrodynamics simula- 
tions. Other studies have confirmed that fragments tend to form near the 



corotation of spiral waves, even when radiative physics is included (Boley and 



Durisen 


2008 


Boley , 


2009 



with a corotation at r, is stable against fragmentation for 
47r2/BHp/9 



5r 
r 



> 



81(sini)2 (Q)m2^ 



(5) 



where m is the number of spiral arms and i is the pitch angle of the spiral. 



which is typically i ~ 10° in gravitationally unstable disks (see Boley and 



Durisen, 2008 Cossins et al. , 2009). Using the DHP2008 definition for /dhp, 



we find that /dhp = HM.'^m/ [ixr), where M. is the Mach number for the 
isothermal shock. Strictly, we are using M. to indicate the density enhancement 
in the spiral arm over the smooth distribution, and this should be kept in mind 
when we quantify our results. If the shock truly is isothermal, then the Toomre 
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wavelength within the overdensity will be 

= 2cl/{GT.M^) = Xt/M\ (6) 

We expect fragmentation to occur only when A'j. becomes equal to 26r. The 
ratio of these quantities is 

Fragmentation should occur at corotation when 

/, ^ 1.8 (g) . (8) 

Now that we have an estimate for the location and radial extent of fragmen- 
tation, we need an expression for the surface density perturbation, relative 
to the axisymmetric density, that will permit fragmentation. We also need to 
know the azimuthal width of the shock, which can be estimated by assuming 
it is similar to the full height {2H) of the disk. For estimating H within a 



spiral arm, we refer to Boley and Durisen (2006), who showed that the scale 



height in the post-shock region of a self-gravitating, isothermal spiral shock 
should be reduced by the factor 



nM)J'-±^] (9) 



due to the increased gas density, where q = external gravity/self-gravity for 
the axisymmetric disk. The self-gravity of a disk at its scale height is well 
approximated by 27rG'S, while the star's gravity by Qcg. Combining these 
terms yields 

q ^ iQ) /2. (10) 
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However, the self-gravity of the unperturbed disk should be included in our 
definition of fg as well. We approximate this effect by setting 

for {Q) < 2 and fg = F{M.) otherwise. Multiple simulations have shown that 



an initial {Q) < 1.4 is required for an isothermal disk to fragment (e.g., Tomley 



etaTI p^94l [Nelson et alTj pS| [Johnson and Gammiej [2003| [Mayer et aL| [20041 



Durisen et al. , 2007). Using {Q) = 1.4 in equation (8) gives fg ~ 2.5, which 
corresponds to ~ 2.7 in equation (11). As {Q) is lowered, the density 
perturbation that is required to induce fragmentation is also lowered. 

Now that we have estimates for the radial extent of the fragmenting region 
(A^), for the width of the spiral {2H, using fg and equation 4), and for the 
surface density perturbation over the axisymmetric disk (EA^^), we can find 
the initial clump mass within the context of a fragmenting spiral arm: 

Mf = 2Xt^. (12) 



Figure 1 shows the initial mass if it were calculated assuming Mt = '7rS(AT'/2)^, 
assuming the Model A temperature profile, and by using our estimate for Mj; 
we assumed (Q) = 1.4. The curve for Model M shows what we expect from 
the model parameters, and the actual clump masses are shown by symbols. CI 
and C2 match our Mf estimate well, and the initial masses are in the gas giant 
regime. C3 also has a mass in the gas giant range, but our Mf calculation is 
an overestimate. This may be due to differences between prompt and delayed 
fragmentation (see DHP2008 for a detailed discussion). CI and C2 form near 
corotation (prompt fragmentation), while C3 appears to form during the col- 
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lision between the wake of CI and an outer arm (delayed fragmentation). In 
contrast to our Mj estimates, indicates that initial masses should be in 
the range of brown dwarfs for r > 80 AU, even for (Q) = 1. This has led to 



recent claims (Kratter et al. , 2009) that, even initially, gas giant-mass clumps 
should be atypical, which is inconsistent with our analysis and our simulation 
data. Figure 2 shows a close-up snapshot of CI before and just after fragmen- 
tation. Fragmentation is clearly confined to the spiral arm. The pre-fragment 
material has an if ^ 1.2 AU and a. 6r ^ 9 AU. For comparison, we estimate 
that H should be ~ 2.2 AU and 6r ^ 8.5 AU using our analysis above. 



4 Initial Angular Momentum 

The radial extent of fragmentation suggests that clumps should have sub- 
stantial angular momentum from the shear in the disk. The specific angular 
momentum of a newly-formed clump can be approximated by 

Jinit ^ l/3(^](r + 5r)(r + 6rf - Q{ry), (13) 

which for 5r <^ r, Jmit ~ l/6{GMsta.r/ry^'^5r, where Sr = A'j./2. This esti- 
mate assumes that the difference between the orbital angular momentum of 
a clump's outermost material and of the centroid of fragmentation (r) goes 
entirely into rotation. The factor of 1/3 comes from assuming we have a rigid 
rod of length 26r for the moment of inertia, based on the shape of the col- 
lapsing region in the left panel of Figure 2. Figure 3 shows the initial angular 
momenta for hypothetical fragments that we expect would form in Models A 
and M, as a function of radius, according to equation (13). Clumps that form 
at r ~ 100 AU should have J ~ few x lO^^cm^s^^, which is roughly consistent 
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with the simulation data. This also strongly suggests that a circumplanetary 
disk can form with a radius of approximately an AU (Fig. 4) after subsequent 
clump contraction. The expected angular momentum radius, i.e., where we 
expect rotation to limit further contraction, can be estimated by 

rj fti . (14) 

36rMf ^ ' 

Using our estimate for the fragmentation mass Mf from section 3, this becomes 

where vk is the Keplerian orbital speed. The sizes of the clumps in the SPH- 
SIM (Table 1) are consistent with this estimate (Fig. 4), even though we 
overestimate J by a factor of about two. This suggests that thermal pressure 
is still an important component in a clump's initial size. Indeed, the initial 
T/|W^| for the clumps is roughly 0.2, where T is the total rotational energy and 
W is its potential energy. These clumps are initially more like rapidly rotat- 
ing spheroids than true central object+disk systems. During the contraction 
phase, the clump may become susceptible to dynamical instabilities, for exam- 



ple, the bar instability (e.g., Durisen et al. , 1986). Such a dynamic event would 



rapidly rearrange the angular momentum distribution of the object, leading 
to rapid outward transfer of angular momentum and overall expansion of the 
clump. Convection could also play a role in redistributing angular momentum, 
but its overall effect is uncertain. Unless some mechanism can transfer angular 
momentum inward from the outer mass shells to the central regions, our es- 
timate should be valid, and shows that the high-J material in the clump will 
be unable to collapse to the size scales of Jupiter. Once dissociative collapse 
is reached, a disk should form. To illustrate this point further, we compare 
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Jupter's rotational angular momentum to a wide-orbit clump's initial J. It 
is unclear whether Jupiter's angular momentum is consistent with the Js for 
planets on wide orbits, but it provides a reference value. For this estimate, we 
assume that Jupiter is a rigidly rotating sphere with a radius ~ 7 x 10^ cm 
and a rotational period of 10 hr. This gives t/jupiter ~ 3 X 10 cm s , which 
is about three orders of magnitude smaller than a clump's initial J. 

We have ignored the possibility that a large fraction of the angular momentum 
given to the clump from shear goes into altering the orbital angular momentum 
of the clump. We expect for this to affect our result by a factor of order 
unity. Likewise, the orientation of the arm just before fragmentation and the 
difference between the assumed moment of inertia (the rigid rod) and the 
actual one will also cause errors of order unity. Nevertheless, the SPHSIM 
data show that our estimate for J is valid to an order of magnitude and 
provides a reasonable upper limit. 

5 Contraction Timescale 

So far, we have defined the mass scales and initial J that we expect for 
fragmentation. Whether these clumps can become bound objects depends, 
in part, on their primary contraction timescale. Once central temperatures 
reach Tc ~ 2000 K, H2 begins to dissociate, and a rapid collapse to sizes of a 
few Jupiter radii will follow because energy goes into dissociation instead of 
thermal support. After the H2 collapse, the clump will be stable against tidal 
disruption owing to its small size. In contrast, a clump that is transported into 
the inner disk before it undergoes dissociative collapse could become tidally 
destroyed. This is illustrated by Figure 5, which shows a snapshot of a clump 
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being tidally disrupted in SPHSIM. Two clumps interact, and C2 from Table 
1 is put on an eccentric orbit. As it approaches periastron, the clump's volume 
becomes much greater than its Hill volume, and it is destroyed. The disruption 
of C2 takes place about 1500 yr after its formation. Is this disruption physically 
motivated? Should the clump have already undergone dissociative collapse? 
In this section, we determine whether the primary contraction timescale (ri) 
is long enough to make clump disruption a real possibility in protoplanetary 
disks. 



Helled (2009, private communication) has graciously shared her data with us 



(from Helled and Bodenheimer, 2009), which show that clumps with masses 
of a few Mj roughly follow an n = 2.3 polytrope. In order to calculate ri for a 
variety of conditions, we have developed a poor-man's gas giant evolution code 
using polytropes. By assuming initial radii and masses as calculated above, 
an initial polytrope solution can be determined. Using the polytrope profile 



and D'Alessio et al. (2001) Rosseland mean gray opacities, a photosphere can 
be calculated, which gives us the luminosity L = AnRl^aT^Q. The time step 
between iterations is At = O.OIU/L, where U is the internal energy of the 
polytrope. For a given step, the total energy of the polytrope is updated ac- 
cording to E = U + W — AtL. The potential energy of the system is then 
calculated as W = 3E-^0g — "^f^^^^, where 7 = 7/5 is the thermodynamic 
adiabatic index, which is not necessarily the same as the structural adiabatic 
index 'jp = 1 + 1/n. The term with AM accounts for mass that is accreted 
over At. Once the new potential energy is determined, a new radius can be 
calculated using R = —3GM'^/{{5 — n)W), and the new internal energy can be 
calculated from the virial, i.e., U = — H^/(3(7 — 1)). With the new mass and 
radius, the density and temperature profiles can be updated. This continues 
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until Tc > 2000 K. For all calculations here, n — 2.3 and the maximum grain 
size assumed in the opacities is 1/xm. Unfortunately, we ignore clump rota- 
tion for these calculations, even though, as discussed above, we expect it to 
influence subsequent clump evolution. To reiterate, this omission is expected 
to have two principal effects: (1) Rotation is extra support, and the contrac- 
tion should be altered, especially when dynamical instabilities set in. (2) A 
nonspherical contraction will likely alter the amount of material in the clump 
that is exposed to high temperatures, but our calculations present an order of 
magnitude estimate. 

We now consider the structure and contraction of hypothetical clumps using 
the method described above. Polytrope Clump 1 (PCI) is followed without 
accretion, and its mass is set to 3Mj with an initial polytrope radius ~ 5 
AU. The photosphere is initially at a radius ~ 3 AU, with a temperature 
of 22 K. The background irradiation could affect the contraction timescale, 
but we do not address that detail here. The central temperature and density 
evolution for PCI is shown in the p — T plane in Figure 6. Symbols indicate 
10^ yr intervals. The contraction time for PCI is 8 x 10^ yr, which is very 
similar to the results of Helled & Bodenheimer. This timescale is about 100 
orbits at 100 AU for a 1.5 Mq star, giving a large window of opportunity for, 
e.g., clump- wave and clump-clump interactions to transport clumps inward 
or outward. However, we also need to address whether accretion will decrease 
the contraction timescale such that disruption becomes unlikely. 

We assume that the clump is accreting at its maximum rate, where M onto 
the clump is limited solely by the rate that mass can be delivered to the 
planet's Hill sphere. For simplicity, consider a cylinder centered on the frag- 
ment, with Hill radius Rh- The flow of material into into this cylinder should 
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be dominated by the shear flow around the fragment, which yields 

r fl (r ) 

M^2 J {r o) Rjj cos i sin idi, (16) 



where Tq is the orbital radius and i is the angle from the perpendicular to 
the radial direction. We have assumed that S, the surface density at ro, is 
constant over perturbations of Rh from tq. Evaluating the integral gives 

M . 2 X 10- yr-' f^)"' (^f f^V'" ^ (17) 

VlOOAUy VMj; \Mq J lOgcm-2'^ ^ 



for clump mass M. We have chosen to normalize the function to E = 10 g cm"^ 
because this surface density corresponds to an unstable disk with Mg^ar = Mq 
and T = 10 K at r = 100 AU. As the clump accretes, the size of its Hill 
sphere grows, allowing it to capture more mass. As it becomes more massive, 
the evolution time toward rapid dissociative collapse decreases. 

Figure 7 shows the ratio of the clump mass just before Tc = 2000 K to the 
initial mass, as a function of the initial fragment mass. Although we expect 
initial clump masses to typically be a few Mj, we show higher masses for 
completeness. Each fragment is assumed to grow at the rate given by equation 
(17), and has an initial polytrope radius of 5 AU. For calculating M, we assume 
a primary mass of 1 Mq, E = 10 g cm~^, and an ro = 100 AU. The accretion 
rate is updated every time step, allowing M to grow as the clump grows. For 
comparison with PCI, a maximally accreting 1 Mj clump (PC2) is shown 
on the p-T plane along with PCI (Fig. 6). This clump grows to about 14 Mj 
before it reaches dissociative collapse. Recall that this mass actually represents 
a rapidly rotating clump, which should form a proto-gas giant/brown dwarf 
+ disk system after dissociative collapse. The timescale to reach collapse is 
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also shown in Figure 7, which indicates that fragments with initial masses less 
than about 6 Mj remain susceptible to tidal disruption for ~ 10 orbits at 100 
AU around a 1.5 M© star. As the initial fragment mass increases, accretion 
becomes less important to the precollapse evolution, but the time available 
for disruption continues to decrease. 

We do not claim that equation (17) represents the rate that a clump must 
accrete material; instead, we use it only as a rough upper limit for the accre- 
tion rate. As discussed above, rotation could affect the contraction timescale 
and lead to dynamical instabilities in the rapidly rotating clump. Additional 
caveats include the following: (1) The material entering the clump's Hill sphere 
may have non-negligible thermal pressure. (2) The energy gained due to the 
accretion of gas may not be efficiently radiated away. (3) Some material en- 
tering the Hill sphere may be unable to shed enough angular momentum to 
become part of the clump. (4) Convection could rearrange the mass and angu- 
lar momentum distribution of the clump, leading to additional instabilities. It 
is unclear how these assumptions will affect out accretion estimate in detail, 
but we expect that these processes will decrease the average mass accretion 
rate of the clump, which will increase the time a clump remains susceptible 
to disruption. We conclude that the disruption of clumps during the earliest 
phases of disk evolution is quite possible. 

6 Consequences of Clump Disruption for Planet Formation 

The work in the previous sections demonstrates that clumps can remain in 
their primary contraction phase long enough after their formation to become 
tidally disrupted. The arguments are based on analytic work as well as simula- 
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tion data. In this section, we discuss three consequences that clump disruption, 
as a general mechanism, can have on disk evolution. Each subsection repre- 
sents a proper study on its own, so our aim here is to present a foundation 
for subsequent work. The first topic, dust processing, is the most strongly 
connected to work presented in this manuscript, while the last topic is much 
more speculative. 



6.1 Dust Processing and Growth 



Figure 8 shows the dust mass in PCI that is above temperatures of 1000 K, 
assuming no mixing or settling and a dust-to-gas ratio of 1/100. Mixing and 
an enhancing the dust-to-gas ratio by enrichment will allow for more solids 



to be processed. Helled and Bodenheimer (2009) show that enrichment by 
planetesimals is unlikely for clumps that form on wide orbits, but this does 
not rule out enrichment at birth. As discussed in section 3, fragmentation is 
most likely to form at the corotation of spiral waves. These spirals will trap 
solids due to the pressure difference between the wave and the surrounding 
disk (Haghighipour & Boss, 2002; Rice et al. 2004), ensuring that the dust 
is enriched where the disk is most likely to fragment. According to Figure 8, 
for roughly 5 x 10^ yr, > IM^ of dust will experience temperatures over 1000 
K for our assumptions. The figure also shows that a clump would need to 
survive for ~ 10^ yr before significant high-temperature, high-pressure dust 
processing is expected to occur. Because C2 only survives for about 1500 yr, 
this particular clump is not expected to process dust thermally, unless it occurs 
in the disruption process itself. However, there could have been substantial 
grain growth during this time, as meter- sized objects can grow in ~ 1000 yr 
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for the conditions in the clump (e,g, WeidenschiUing, 2000). Once the clump is 
disrupted, this material will be distributed into the disk. Because a clump on 
an eccentric orbit can experience multiple phases of mass loss due to changes 
in its Hill sphere, complete disruption is not necessarily required to liberate 
processed dust. 

For clumps that do experience high temperatures, their disruption, whether 
partial or entire, may be a mechanism for processing Calcium-Aluminium- 
Inclusions (CAIs) or CAI-like particles. CAIs are composed of refractory min- 
erals that are stable at temperatures above 1400 K (e.g., Scott 2007). The 
advantage of this formation mechanism is that it can occur as soon as the 
disk begins to form; CAIs are among the oldest objects in the Solar System. It 
is also a fallacy to assert that all of this material will be lost due to accretion 
onto the star. GIs are known to transport material over large radii inward 



and outward, with significant mixing and stirring (e.g.. Boss, 2004; Boley and 



Durisen, 2006). The retention of such material should be addressed, which we 
leave to future work, but we do not expect for all of this processed material 
to be lost. 



6.2 Core Formation 



Clump disruption also may be capable of producing Earth-mass cores early in 
a disk's evolution. These cores could provide a substantial head start for the 
core-accretion mechanism, and are not subject to the meter-barrier problem. 



To illustrate this possibility, we refer to Helled and Schubert (2008), who found 
that cores ~ IM0 can form in contracting proto-gas giants. In their models, 
core formation was halted when the core temperature reached 1300 K, which 
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begins silicate sublimation for the small grain sizes that they considered. For 
clump masses > 5Mj, they found that core formation was inhibited due to the 
fast evolution and high temperatures in the core. We note that larger grains 
and highly refractory materials can survive beyond 1300 K (see, e.g., Fig. 2 in 



Scott, 2007), increasing the time (mass range) available for core formation. 



Any core that forms from a failed gas giant will not necessarily be lost due 
to type 1 migration, although it is not expected to stay at the disruption 



radius. Laughlin et al. (2004) showed that the magnetorotational instability 
can lead to density fluctuations that will put cores on a random walk through 
the disk. GIs should cause a similar behavior, but this needs to be explored 
in future work. The cores that form in failed gas giants could become the first 
planetoids formed in planetary systems, and become seeds for the growth of 
gas and ice giants, or left as dwarf planets in an outer disk. Earth-size objects 
in the Kuiper belt would be consistent with clump disruption. 



6.3 FU On 



A massive clump will remain susceptible to tidal disruption for 10^ to 10^ 
yr. During this time, transport of a clump into the innermost regions of the 
disk (r < 1 AU) may become possible through clump scattering and/or disk 
torques. This topic warrants further study, but we only emphasize this as a 
real possibility during the early phases of disk evolution while heavy infall is 



present on the disk (jVorobyov and Basu| 2006 Boley| 2009|). For a ~ O.3M0 



star, the Hill radius of a 3 Mj object is ~ 0.3 AU at a disk radius of 2 AU. For 
our simple contraction model, we find that the radius of PCI is comparable 
to this size just as it hits ~ 2000 K. Clump disruption inside a few AU 
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can suddenly supply the inner disk with mass, which may lead to a thermal 



instability (Bell and Lin, 1994), and consequently, an FU Orionis outburst. 
Note that we do not require the clump to make it all the way to the star. For 
a very eccentric scattering event, such that the clump arrives at perihelion at 
nearly the escape speed, a Jupiter-mass fragment around a 0.5 Mq primary 
will cause a maximum change in the primary's radial velocity of about 60 
m/s. Unfortunately, the radial velocity drifts of, e.g., FU Ori have only been 
constrained to 300 m/s (Petrov & Herbig 2008). Larger clump/primary ratios 
and smaller perihelion passages may allow for an observable radial drift in 
other systems. 



7 Conclusions 



We have demonstrated that initial clump masses are expected to be in the 
gas giant regime. Our analytic mass estimates are consistent with hydrody- 
namics simulations with radiative cooling. These clumps will have significant 
angular momentum, suggesting that disks ~ AU in radius should form after 
dissociative collapse; the initial clump mass represents what will eventually 
become a core-disk system. Even when mass is assumed to be accreting at its 
maximum rate, at least with our analytic estimates, 1 Mj clumps are only 
expected to grow to masses ~ 10 Mj before H2 dissociation causes rapid col- 



lapse, consistent with known masses for the planets in HR 8799 (Marois et al. 



2008). Further growth of the system may be regulated by the circumplan- 



etary/brown dwarf disk due to the angular momentum of the new material 
entering the Hill sphere. Clump ti contraction timescales, even for several Mj, 
will be 10^ to 10^ yr, giving sufficient time for clumps to be transported into 
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the inner disk and to be tidally destroyed. Clumps will have very different en- 
vironments from the typical conditions in the outer disk, and they represent 
factories for processing dust and building large solid bodies. Clump disruption 
therefore represents a mechanism for processing dust, modifying grain growth, 
and building large, possibly Earth-mass, objects during the first stages of disk 
formation and evolution. 
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Clump Mass (Mj) -Rgphere (AU) a (AU) J (cm^s""^) Formation (yr) 



CI 



3.3 



2.9 



80 



2.3(18) 



300 



C2 



1.7 



2.3 



65 



1.1(18) 



750 



C3 



0.87 



1.8 



85 



3.4(17) 



970 



Table 1 



Initial clump parameters from the simulation SPHSIM. The value -Rgphere is the 
radius that a sphere would have for the same volume as the clump. The formation 
radius is indicated by a. The masses and Js for CI and C2 are consistent with our 
estimates in sections 3 and 4. C3 seems to form as a result of an interaction between 

an existing clump's spiral wake and a wave, and is consistent with delayed fragmen- 
tation. The last column indicates the formation time after the start of SPHSIM. 
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Fig. 1. Mass estimates calculated by assuming the Toomre mass Mt = 7rS(AT/2)^ 
and by calculating Mf as described in the text. The curve associated with Model 
M represents the fragmentation mass that one would expect in the simulation. The 
curve labeled "T(K) Model A" gives the temperature profile for Model A. The 
simulation data (crosses) are in good agreement with the Mf estimates. 
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Fig. 2. The particle distribution for CI just before and after fragmentation. Green 
represents particles that are included in the initial mass given in Table 1. 
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Fig. 3. Specific angular momenta for Models A and M. Model M shows what we 
expect to see in SPHSIM. The simulation data are indicated by crosses. 
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Fig. 4. Angular momentum radius rj for or Models A and M. Model M shows 

what we expect to see in SPHSIM. The simulation data for i^sphcrc ^ire indicated by 
crosses. Because the Js for the clumps are lower than expected by a factor of order 
unity, the correspondence between -Rgphere fJ be due to thermal support. 
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AU 



Fig. 5. Surface density snapshot from the simulation SPHSIM. The innermost clump 
(now quite elongated at x,yRi20,5 AU) is being tidally destroyed, and does not 
survive for more than 1/4 more of an orbit. CI, near the top of the snapshot, is 
becoming bar unstable. 
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Fig. 6. The p — T plane for the contracting polytrope models. Symbols indicate 
roughly 10^ yr intervals. Once the central temperature reaches 2000 K, rapid collapse 
should occur. PCI is a 3 Mj clump contracting in isolation. PC2 starts at 1 Mj 
and accretes mass as fast as gas can flow into the clump's Hill sphere. 
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Fig. 7. The ratio of M2000 to the initial fragment mass Mi, where M2000 is the 
mass of the clump when the central temperature reaches 2000 K, leading to rapid 
collapse. The contraction time in kyr is shown with the dashed line, and is on the 
same scale as the mass ratio curve. Even with accretion, clumps with initial masses 
below ~ 6 Mj remain susceptible to disruption for 10^ yr. For masses greater than 
roughly 8 Mj, the clump collapses before the mass can be doubled. By Mj = 20 
Mj, accretion is marginalized. 
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Fig. 8. Dust mass in regions with T > 1000 K. We have assumed no setthng or 
mixing and a dust-to-gas ratio of 1/100. 
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